Hamiltonian two-body system in 
special relativity 

<n : 

q3 1 Philippe Droz- Vincent 

LUTH 

m ■ 

(N ; Meudon * 

S 1 : 

<— i 

^3 Abstract 

We consider an isolated system made of two pointlike bodies inter- 
'— — acting at a distance in the nonradiative approximation. Our framework 
\ is the covariant and a priori Hamiltonian formalism of "predictive rel- 

■ ativistic mechanics", founded on the equal-time condition. The center 

. of mass is rather a center of energy. Individual energies are separately 

^-j. | conserved and the meaning of their positivity is discussed in terms of 

' world-lines. Several results derived decades ago under restrictive as- 

\Q ■ sumptions are extended to the general case. Relative motion has a 

structure similar to that of a nonrelativistic one-body motion in a sta- 
tionnary external potential, but its evolution parameter is generally not 
r • ■ a linear function of the center-of-mass time, unless the relative motion is 

• i-H , circular (in this latter case the motion is periodic in the center-of-mass 

^ | time). Finally the case of an extreme mass ratio is investigated. When 

this ratio tends to zero the heavy body coincides with the center of mass 
provided that a certain first integral, related to the binding energy, is 
not too large. 



1 Introduction 

Classical relativistic dynamics of pointlike bodies has a long story; with- 
out claiming to be exhaustive let us mention the Wheeler-Feynman 
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(WF) electrodynamics [JJ based upon the Fokker action [2], the three 
forms of dynamics (front form, point form, instant form) advocated by 
Dirac [3j , and after the discovery of a famous No-Interaction theorem [3] , 
various efforts made in order to circumvent it; for instance Predictive 
Mechanics [5] [6], the Singular Lagrangian method [8] and Constraint 
Dynamics [9j. In the last decade there were a few papers along the lines 
of WF [10] and also the work carried out by Lusanna et al. [TTJ in order 
to give a covariant status to the instant form. 

Beside several n-body generalizations, most progresses have been de- 
voted on the two-body problem, which is our present subject. 

A first point was the possibility of actually having (unlike WF elec- 
trodynamics) second-order differential equations describing the motion 
of a finite number of degrees of freedom interacting at a distance. Ac- 
cording to this view, the field that carries interaction is supposed to be 
eliminated, the space of initial data has a finite number of dimensions 
and a point in this space uniquely determines the subsequent motion [5] 

A second point was about a Hamiltonian formalism. We have in mind 
conservative mechanics: radiative corrections are neglected and we fo- 
cus on isolated systems, characterized by a finite number of degrees of 
freedom and by Poincare invariance. Under these conditions the no- 
interaction theorem [3] excludes the possibility of demanding that the 
physical positions be canonical throughout phase space (as was always 
done in classical mechanics). Relaxing this requirement implies an ar- 
bitrariness which can be removed by imposing a relationship between 
physical positions and canonical variables on some submanifold S. 
In the a priori Hamiltonian approach of predictive mechanics [7] , man- 
ifest covariance is realized with help of considering degrees of freedom 
that are geometrically natural but redundant (if compared to the clas- 
sical situation). Positions and momenta are four-component objects, 
phase space has sixteen independent dimensions and we employ a mul- 
titime formalism (in a different spirit Todorov proposed to focus on the 
physical degrees of freedom only, in another formulation of dynamics 
which is covariant as well but includes first-class constraints [9]). 

After discussing other possibilities, we put forward a natural pre- 
scription: physical and canonical positions must coincide when both 
physical positions are simultaneous with respect to the rest frame of 
the system [7]. In this approach a submanifold S is obtained by se- 
lecting the configurations where the physical coordinate times x\,X2 
are equal in the center-of-mass frame. An advantage of this equal-time 
prescription is to permit the contact with the constraint formalism, as 
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shown in detail by L. Lusanna [12j . 

The two-body problem immediately suggests these important issues: 
center of mass and relative motion. In this article we intend to study 
their properties in the context of relativistic dynamics, aiming at a 
method in order to simplify the determination of the world lines. 

A few exploratory results that we derived in the past [33] |14j have 
remained fragmentary, as most of them have been obtained with help of 
restrictive assumptions concerning the shape of the interaction poten- 
tial. The present work is free of such limitation, aiming at the possibility 
of dealing with realistic interactions; such interactions cannot offer the 
simplicity of the academic models we considered long time ago: see for 
instance the unavoidable P 2 - dependence of an electromagnetic two- 
body potential proposed by Jallouli and Sazdjian [20], in the framework 
of relativistic quantum mechanics. 

Before we focus on its equal-time version, we shall sketch the main 
lines of the a priori Hamiltonian approach in general. 
In the next section we present the isolated two-body systems and intro- 
duce the center of mass. 

In the rest of the paper we consider unipotential models (for these sys- 
tems the individual energies are separately conserved). 
In Section 3 we study the evolution of the canonical variables, irre- 
spective of their relationship to physical positions. Later, in Section 4, 
we focus on the equal-time prescription and discuss a strategy for the 
determination of world-lines; the importance of circular orbits is em- 
phasized. In 5 we consider the case of an extreme mass ratio; a toy 
model is presented in Section 6, and Section 7 is devoted to concluding 
remarks. 

The velocity of light is taken as unity, except in the Appendix. When 
no confusion is possible, tensor indices are omitted and the contraction 
dot is employed also for tensors, for instance JP stands for the vector 

1.1 The a priori Hamiltonian formalism 

The canonical coordinates are q^^Pbp i n the sixteen-dimensional phase 
space; q\,q2 are points in Minkowski space M. and p\,P2 are four- 
vectors. 

The symplectic form dqf A dp\ a + dq% A dp2 a corresponds to the Poisson 
brackets 

{<£,Pbf}} = 6abtf (1) 
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Note that the particle labels a, b, c are not summed over when repeated. 
The Hamiltonian equations of motion 

^ = {q a ,H b } ^ = {pa , Hb } (2) 

involve two Hamiltonian generators Hi , H 2 submitted to the predictivity 
condition 

{H 1 ,H 2 } = 

The two evolution parameters t\ , t 2 are suitable generalizations of the 
proper times (normalized to the masses). Our notation is choosen such 
that the generators can be identified as the half-squared masses; this 
allows to consider the masses as first integrals, redeeming the fact that 
phase space has redundant degrees of freedom. 

For instance, in the trivial case of two free particles one is left with 

2H a = p 2 a , the physical positions reduce to x a = q a and the evolution 

parameters are just r a = s a /m a , where si, s 2 are the arc lenghts. 

For interacting particles the Hamiltonians involve additional terms V\ , V 2 referred 

to as "potentials". In general qi,q 2 differ from the physical positions, 

in fact x\,x 2 are determined by the partial differential equations 

{x a ,H b } = for a^b (3) 

and some reasonable initial data. Solving these equations provides a 
correspondance between physical and canonical coordinates, say 

x, v < — > q,p 

which, inserted into the solutions of system ([2]), ensures that x\ is a 
function of T\ only and x 2 is a function of r 2 only, in other words ([3]) 
ensures the existance of one-dimensional world-lines [7]. This point is 
easily checked by introducing on phase space the Hamiltonian vector 
fields X\ , X 2 defined by 

X a f = ^- = {f,H a } (4) 

for all phase-space function /. The predictivity condition corresponds 
to the vanishing of the Lie bracket [Xi,.^]. Although the analytic 
shape of the Hamiltonians yields some information about first integrals, 
symmetries, etc, the physical meaning is fixed only once a solution to 
(|3|) has been specified. At this stage phase space is identified with the 
bundle product T(M) x T(M). 

Notation. (H a ) iree = \p\ , they generate (X a ) ime . 
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When fixing numerical values \m\ and \m\ to the Hamiltonian gener- 
ators, it is convenient to set 

^=\im\+m 2 2 ), y =\ijn\-m\) (5) 

which amounts to 

mim2 = J u 2 — v 2 , {mi + "I2) 2 = 2// + 2y a 2 — v 2 (6) 



2 Isolated two-body systems 

The Lie algebra of the Poincare group is spanned by P a , J^ u ■ We sepa- 
rate external from internal variables by setting 

p = pi+P2 Q = \{qi + q2) (7) 

y=^(Pi-P2) z = q 1 -q 2 (8) 

Remark in general Q is not the center of mass El- 
Other definitions of Q and y (as conjugate to P and z respectively) were 
possible, but the advantage of ours is that it does not require to a priori 
fixthe numerical values of the constituent masses. 

The standard Poisson brackets that do not vanish can be re-arranged 
as follows 

{Q<*,Pp} = {z<*,yp} = 6% (9) 

and we can write 

J = QAP + zAy = q 1 Ap 1 + q 2 /\p2 

Naturally Hi, H 2 are supposed to be Poincare invariant (vanishing Pois- 
son brackets with P and J) and the initial conditions for solving ([3]) must 
be invariant as well. 

From we can derive several useful formulas; let us list some of 
them. Of course {Q^,P 2 } = 2P^ and 

so if we define the projector orthogonal to P 

n-,-^ do) 



^ven for free particles it is center of mass only when the masses are equal. 
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we get 

{Q^z 2 } = -2z^ (11) 
with this notation that the tilde means an application of II, say 

r = n a/3 fa 

the r.h.s. of (|lip is orthogonal to P; it follows that 

{Q ■ P, z 2 } = (12) 

On the other hand we compute 

P (5?) P 

{Q ■ p, -^-} = (13) 

in other words Q ■ P has a vanishing Poisson bracket with the projector 
IT. Moreover we easily check that 

{Q ■ P, {JL ^} = (14) 

It is obvious that {Q-P, y 2 } vanishes and we can apply the above formula 
to the identity 



~2 2 

y =y 



(y ■ P? 



P 2 

which yields 

{Q-P,y 2 } = (15) 

For the relative variables note that the spatial piece of the one has a 
vanishing bracket with the time piece of the other (its conjugate): 

{z, y-P} = {y, z-P} = (16) 

Similarly it stems from (|13p that 

{Q-P,z a } = {Q-P,y a } = (17) 

Of course we already know ten first integrals, namely P a , J^ v . We 
obviously have that 

-752- = ~pT y ~ Q ~ -752- z ( 18 ) 
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2.1 Center of Mass 

The possibility to define a center of mass, using linear and angular 
momenta as basic ingredients, has been known long time ago [16J. 
In a previous work |15j we proposed this canonical definition for the 
components of the center of mass 

S = ^ + (^#)P (19) 

or equivalently 

z = Q + ( y -^r)z-( z -j£)y (20) 

This can be transformed again, if we notice that 
P 2 /2 db y ■ P = P - pi, (resp. P ■ P2), we get 

„ P-Pl P-P2 P-Z . . 

a = -^-9i + -^92 - -^y (21) 
Formula (|20p entails 

E-P = Q-P (22) 

In (|19p the only quantity which is not a constant of the motion is 
(P • Q/P 2 ) , it has the same dimension as t\,T2. In contrast setting H. 

T= p -Q 
\p\ 

we give to T the dimension of lenght (and time since c = 1). One easily 
computes 

whence we derive the relations 

{E a , P p ] = 5<$ ^} = P a (23) 

Owing to the constancy of J^ u and P a , we can fix these quantities, say 
in particular 

Pa = k a , timelike vector 

define k a k a = M 2 

Then we see that formula (|19j) defines the coordinates of a point which 
moves on a straight line when both t\,t% run freely and independently 
from —00 to +00, and the direction of this line is given by k a . 



2 definition of T has a different dimension in [TS] , but the same as here in [14] 
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P ■ Q k a 

Equation (fT9l) becomes E a = Const. + ( — — — ) — . 

— v M M 

We can write —j^- = — and consider respectively T and T/M as the 

proper time and the evolution parameter of the center of mass. 

Similarly equations ([231 ) supplemented by the trivial observation 

that {P, 2^ >2 } vanishes, can be viewed as canonical equations of motion 

for E, generated by the one-body Hamiltonian ^P 2 . 

Note however that in (|19j) the components of E do not commute 
among themselves; a similar situation was already encontered by Pryce [T6] . 

From now on we shall focus on Unipotential Models characterized 
by the same interaction term for both particles, say V\ = V 2 = V . 
Such models permit to write down explicit forms of the Hamiltonians. 
Moreover it seems that they are still general enough for dealing with 
most realistic interactions. 



3 Evolution of the canonical variables 

For the moment we postpone the physical interpretation; we are inter- 
ested in the evolution of the canonical variables in terms of the param- 
eters t\,t 2 , more generally we are concerned with statements that hold 
true regardless to the prescription used for solving the position equations 
([3]). In this section, the only thing we assume about this prescription is 
Poincare invariance. 

At this stage the separation of external/internal variables is just 
formal and convenient for easy calculations. Because of Poincare invari- 
ance and predictivity, V can be only a function of the five independent 
scalars [13] 

P 2 ,I 2 ,y 2 ,z-y,yP (24) 

But in view of ()14j) it is more convenient for calculations to re-arrange 
them as 

P ,z ,y ,z-y, 2 (25) 

Unless otherwise specified we consider V as a function of these five 
arguments. Soon it was observed [13] that 

(X 1 -X 2 )z = {X 1 -X 2 )y = 

implying this 

Proposition 1 In the motion, z and y will depend only on t\ + t 2 . 
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In contrast, z ■ P and Q ■ P will depend on both evolution parameters. 



For any unipotential model, the individual energies ^ , ^J^ 2 

are separately conserved. Moreover the translation invariance of V im- 
plies conservation of P 2 , hence this obvious first integral 

N = H 1 + H 2 - ^- p2 H ^ (26) 

By elementary manipulations we find 

N = y 2 + 2V (27) 

This important function defined on phase space is intimately related 
with the properties of relative motion; this can be intuitively seen as 
follows: 

fixing numerical values to Hi, Hi, P a (in particular P ■ P = k ■ k = M 2 ) 
results in a numerical value for N, let it be 

< N >= -A 

The quantum mechanical analog of this quantity appeared, denoted as 
b 2 , in the work of Todorov |19j . 

Employing (for instance) the reduced mass of Galilean mechanics, say 

mi7B2 

m 



mi + rri2 

we can start from (|27p and check that A/2mo is the leading term in 
the post-Galilean development of the quantity M — (mi + TO2) which is 
usually considered as binding energy for the bound states of relativistic 
quantum mechanic (see Appendix I). 

In addition we shall see later on that generates the evolution of the 
spatial relative canonical variables according to a one-parameter Hamil- 
tonian scheme reminiscent of the nonrelativistic one-body mechanics, 
see equations ([35]) below. 

From (|26p and with the notation ([S]) we can write 

M 2 v 2 

A = — + W 2 ~» (28) 

This important relation between A and M can be solved for M 2 . We 
first write 

M 4 - 4(/i + A)M 2 + Au 2 = (29) 
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From (|28p it is already clear that A + \x > 0. But M 2 must be real 
and positive, so the possible values of A are further restricted by the 
condition 

\H + A\>\u\ (30) 
which ensures that M 2 is real; under this condition we can write 



M 2 = 2(/j + A)±2 JGu + A) 2 -!/ 2 (31) 



Moreover we must have at least one root of (J29J) positive. Since their 
product is non-negative these roots cannot have opposite signs, thus 
it will be sufficient to ensure that their sum is positive. Hence the 
condition 

/i + A>0 (32) 

which ensures M 2 > 0. We can encompass both (j30|) and ([32]) by 
writting 

H + A>\v\ (33) 



Fortunately, the sign ambiguity in (|3ip can be removed, with help of 
the individual positive-energy condition (remind that P ■ p a /M are the 
individual energies) that we assume henceforth 

P ■ pi > 0, P ■ p 2 > 

M 2 

Indeed, the numerical values of P ■ p\ and P ■ p2 are given by — — ± v. 
Requiring that both are strictly positive amounts to the condition 

M 2 > 2\u\ (34) 

Let M' 2 , M" 2 be the roots of equation ()29|) . If both roots were to satisfy 
this inequality, it would contradict the equality M' 2 M" 2 = Av 2 implied 
by the last term in (|29p . It follows that 



Proposition 2 Under the condition of positive individual energies, only 
one root of A29\) is admissible 



( Remark For strictly equal masses, only the plus sign may be taken 
in (|3"Tjh Indeed this case means v = 0, one root of equation (|2Uj) is 
obviously zero which must be rejected, and the other root is M 2 = 
4(^i + A) obtained from (|3~T1) by choosing the plus sign). 

We can assume that m\ < m,2 without loss of generality. Thus v < 0, 
and (|33"j) becomes 

fi + v + A > (35) 
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or equivalently 

m\ + A > (36) 
We see that either A > or it satisfies — m\ < A < . In other words, 

Proposition 3 Under the assumption of positive individual energies, 
either A > or |A| <m\. 

1 2 

If we impose the individual energy conditions we additionally get — M > 
v (trivial) and also -^M 2 > —v, in other words 

M 2 > 771-2 — m i (37) 

Thus, as soon as mi ^ 771,2, the collective mass of the whole system 
(lenght of linear momentum) cannot be arbitrarily small. 

Let us check that taking the plus sign in (1311) always yields an admissible 
root; eqn ([35]) implies that \x + A > — u, thus looking at (f3T]) we can 
write 

— > -v + \J(n + A) 2 - v 2 > -v 
implying the individual energy condition (JMJ). [j 

In view of the Proposition 2 above, taking the minus sign in ([31 p is 
excluded and we can write 

M 2 = 2(/i + A) + 2 + A) 2 - v 2 (38) 

Note that M reduces to m\ + m<i in the nonrelativistic limit, see Ap- 
pendix I. 

We now investigate the evolution of the dynamical variables. Let us 
first analyze the evolution of z and y. 

(X! + X 2 ) z a = {X l + X 2 ) irce T + 2{z a , V} (39) 

(X, + X 2 ) y a = (X, + X 2 ) ime y a + 2{y a , V} (40) 
where V is a function of the quantities listed in (|25l) . 

Proposition 4 T/ie Poisson brackets {z a ,V} and {y a ,V} are combi- 
nations of z a , y a with coefficients that are functions of the five scalars 
listed in [25\) . 
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Proof. 

Consider first z a . Obviously {z a , P 2 } = {z a ,z 2 } = 0, and we also have 
that {T,yP} = 0. 
Then we compute 

{I a ,y 2 } = 2y a 
{z a ,z • y} = z a 

hence 

dV dV 
{ »V- ) = 2 _ r + _» (41, 

Since V is a function of the scalars (|25p only, its partial derivatives 
involved in (|4ip obviously share this property. 
Then consider y a . Obviously 

{r,P 2 } = {y a ,y 2 } = {y a , ijL ^} = o 

Then we compute 

{y a ,z 2 } = -2z a 

{y a ,z-y} = -y a 

hence 

dV dV 

^ = -^-Wl) r (42) 

The partial derivatives involved in this formula are functions of the 
scalars 



In view of Prop. 1 it is convenient to set 

A = n + r 2 (43) 
the equations of motion for z, y are 

§-«^ + V) (44) 

t = iv\? + V} (45) 

where the brackets can be computed as functions of z^,y u , and of the 
first integrals P 2 , y-P ■ Once P 2 and y-P have been fixed, the evolution 
of the spatial internal variables is given by a system of six first-order dif- 
ferential equations, to solve for six unknown functions; this problem has 
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the structure of a nonrelativistic problem for one body in three dimen- 
sions. The four-vectors z and y remain within the 2-plane orthogonal 
to k and to the (conserved) Pauli-Lubanski vector [13]. Some solution 



of the system (j44]) (|45p defines the evolution of the spatial relative canon- 
ical variables z, y. Interpretation in terms of world lines, relative posi- 
tions and relative orbit will be given in the next section. 
The collective evolution parameter A plays the role of the Newtonian 
time in the analogous one-body system. But in general A is not the time 
of any inertial observer. Therefore, in order to evaluate the schedule of 
the relative motion, we should express A in function of T and insert the 
outcome into £. 

A complete knowledge of the motion also requires that we determine 
the evolution of z ■ P and Q ■ P in terms of t\,t 2 . 
By sum and difference 



We know that V depends only on the five scalars (|25p and we observe 
that {z ■ P,y a } = 0, implying that z ■ P has a vanishing Poisson bracket 
with all scalars ([25]) except (y ■ P) 2 /P 2 ). We find {z ■ P, y ■ P} = P 2 
hence finally 

{z ■ P, V} only depends on the five scalars ([25|h 
After integrating the system ([4"1"|) (ji"5j) let us set 



z = ((\,P 2 ,y-P) 



y = V (X,P 2 ,y-P) 



{X 1 +X 2 ) z-P = 2y-P + 2{z ■ P, V} 



(46) 



{X l -X 2 )z-P = P 2 



(47) 



{z-P,V} = G(\,M 2 ,v) 



(48) 



We obtain 



{Xt + X 2 )z -P = 2v + 2G(A, M 2 , v) 
(Xi - X 2 )z ■ P = M 2 



Since 




(X 1 -X 2 ) = A_ JL = 2 

OT\ OT 2 



d(n - t 2 ) 







we finally have 




(49) 
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Observing that 

M 2 

—y- ± v = P • pi (resp.P • p 2 ) 
we may write equivalently 

z-P = (P- Pl ) Tl - (P ■ p 2 )r 2 + J GdX + const. (50) 

which reduces to eq. (3.6) of [13J when G vanishes. 

Similarly in view of (|14p and (|17p we can simply write (with d according 

to d25]) ) 

r)V 

where {Q ■ P, P 2 } = 2P 2 therefore 

f)V 

{Q-P,V} = 2P 2 ^ (51) 

which only depends on the five scalars (|25p. 
After integration of the system (|44l) (l45]l let us set 

{Q-P,V} =F(\,M 2 ,v) (52) 

in other words we perform the substitution 

F = subs. (z = C, y = rj, P 2 = M 2 , y ■ P = v\ {Q • P,V} ) (53) 



Now we can write 

(X 1 -X 2 )Q-P = {Q-P,yP} = yP (54) 

(X 1 + X 2 )Q ■P = (X 1 + X 2 ) bee + 2{Q • P, V} (55) 

(X 1 + X 2 )Q .p=±P 2 + 2{Q • P, V} (56) 

Straightforward integration yields 

v M 2 f 

Q-P = -(t 1 -t 2 ) + —A + Fd\ + const. (57) 

We see that modulo the solving of |^[ ) |^5| j the evolution of Q ■ P in 
terms of t\,t 2 will be given by a quadrature. In the very special case 
where F = (with our present notation) the above formula reduces to 
eq (3.8) of [13]. But most realistic potentials actually depend on P 2 
which implies that F differs from zero. 
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In contrast G vanishes in several cases of interest, for instance 
G = provided V depends only on the dynamical variables z 2 , P 2 , L 2 . 
This statement stems from (|16p . 

To summarize: After integrating (|44 p(|45p we got z ■ P and Q ■ P 
as functions of t\ , t% . Since we have the first integrals P a = k a and 
y ■ P = v the only remaining dynamical variables to be determined are 
QP. According to ([Tgj) 

^ _ z-P _ y ■ P _ PTk 

where the last term is also a first integral, thus everything in the right- 
hand side is already detrmined. This observation achieves to determine 
the evolution of qi,q2,Pi,P2, say 

q a = <j>a(n,T 2 ), Pb = ipb(Ti,T 2 ), (58) 

These functions define the two-dimensional "orbits" El of the evolution 
group in phase space. 

4 World lines in Unipotential Models 

The solutions of system (|2|) can be interpreted in terms of world lines 
provided we ultimately introduce the physical coordonates 
functions of the canonical coordonates. Our Cauchy surface for solving 
the position equations ([3]) is (£) defined by P ■ z = , and our initial 
condition 

x a~qa = ® on the- surface (S) (59) 

can be formulated also as 

x a a =q a a +0{P-z) (60) 

where 0(P ■ z) symbolically represents any expression which vanishes 
with P ■ z. In principle formula (j58j) must be inserted into the solutions 
of ([3]) and yields the worldlines in terms of the individual evolution 
parameters, say x\{t\), X2(t2) ■ 

The above prescription offers several advantages. First of all, setting 



3 To avoid confusion with trajectories in space, we put the word orbit between quotation 
marks when it is meant in the group-theoretical sense. 
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condition (|59p implies that also P ■ r vanishes on (£). In the rest frame 
we can write x\ = x 2 = T, thus finally the manifold (£) can be called 
the Equal-Time Surface . 

At equal times the radius- vector r a moves on a curve that we may call 
the relative orbit. As observed long time ago [13J, this curve lies on the 
2-plane mentioned in the previous section (orbital plane). 

This version of the Hamiltonian formalism clarifies the formal defi- 
nition (|20p written for the center of mass. 
Indeed ([20]) is equivalent to (f2T|) . say 

E= {P ' Pl) gi + ( p '**) ** + o iP . z) 

In terms of the individual energies M a = (P ■ p a )/\P\, we have 
„ M x q x + M 2 q 2 



Mi + M 2 



+ 0{P-z) (61) 



Now at equal times P ■ z vanishes; fixing P a = k a we can replace q a by 
x a and P 2 by M 2 , so we are left with 



Mixi + M 2 x 2 
Mi + M 2 



(62) 



Notice that Mi + M 2 = M and remember that the individual ener- 
gies reduce to the masses in the nonrelativistic limit (see Appendix II). 
In view of these remarks, formula (|62p is more intuitive and significant 
when P ■ pi and P ■ p 2 are both positive: the analogy with the New- 
tonian definition of center of mass becomes obvious, which legitimates 
the positive-energy condition. At this stage it is clear that definition 
(|20p agrees with the one proposed by Pryce [16] ; similarly, formula (j62[) 
agrees with the notion of center of energy according to Fischbach et al 



4.1 Rest-Frame description 

In practice, instead of trying to describe the motion in terms of the in- 
dependent parameters t\,t\ we have better to fix the linear momentum 
k a , so defining a slicing of spacetime by the three-planes orthogonal to 
k a . These three-planes intersect both world-lines, which provides the 
rest-frame description of dynamics, as follows: 

among all possible couples x\,xi the slicing selects the equal-time 
configurations, characterized by k-r = 0. We just have to determine the 
sequence of these configurations, that is a one-parameter set. Picking 
up the equal-time configurations obviously induces a relation between 
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7~i and T2 , by cancellation of P ■ z in formula (|49p . One is left with a 
(possibly nonlinear) expression of t\ — ti as a function of A. 

After imposing (|59p . formula (|49p permits us to express everything in 
terms of A only. Formula (|20p permits to write 

Q2 = ^-(^ + \) z + 0(P-z) (64) 

But z = z + 0(P ■ z). Taking (|59H into account (which implies z = 
f + 0{P ■ z)), we put P ■ z equal to zero in (|63p (|64p and obtain the 
equal-time description of the motion 

Xl =E-{^- l -)a\M\v) (65) 

z 2 = H-(^ + i)C(A,MV) (66) 

These formulas yield a representation of both world lines in terms of the 
same parameter A. But for the sake of a better understanding of the 
motion it is interesting to express A as a function of the center-of-mass 
proper time T. 

Cancelling P ■ z in (I49p . our equal-time prescription implies 

^M 2 (ti - r 2 ) = -uX - J GdX + const. 

Inserting into f)5T[) yields the common value of H° and Q° in the center- 
of-mass frame, say 

T = A (f " W ] - W I GdX + 17i I FdX + const - (67) 

In principle this equation must be solved for A. Let us stress that only 
in the very special case where G and F are constant, the time of the 
center of mass is for all orbits a linear function of the parameter A . 
For instance, this situation is realized when both dV/dP 2 and dV/d(y ■ 
P) identically vanish, implying G = F = 0. This situation will be 
referred to as the Academic Case. 

Otherwise, in the most general case A and T are related in a nonlinear 
way, but for exceptional orbits. 

Note that for a physically admissible solution to (I44h (J45]), T should 
monotonously increase as a function of A. In the academic case this 
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is automatically ensured by the condition (|34p of positive individual 
energies. When the interaction potential V is more complicated we 
must demand dT/dX > 0, whereas we can write 

dT _ M _ v*_ _ uG_ F_ 

dX ~ 4 M 3 M 3 + M ( ' 

For example assume for a moment that G = 0, we are left with a simpler 
condition. If F is positive no problem; otherwise the positive-energy 
condition must be replaced by a more restrictive and model-dependent 
condition (e.g. see the toy model of Section 6 ). 

More generally, the discussion remains easy when F and G are bounded; 
as we shall see in the following section, this circumstance arises in case 
of circular motion. 



4.2 Circular Motion 

Circular motion is characterized by the constancy of z 2 , which implies 
that at equal times r 2 also is constant. We can check that 

Proposition 5 On any circular orbit the functions G, F and the five 
dynamical variables [25\) are constant. 

Proof. 

The interaction potential is a function of the five dynamical variables 
(|24p We have four remarkable constants of the motion P 2 ,N,y ■ P and 
the square of angular momentum 

L 2 = z 2 y 2 -(z-y) 2 (69) 

when fixed, they respectively take on the following numerical values 

M 2 ,-A,u,l 2 

In the set (I24p we can replace z 2 , y 2 , z ■ y by the equivalent set of scalars 
z 2 ,^ 2 ,^. So let 

V = f(P 2 ,z 2 ,y 2 ,L 2 ,yP) 
Since N = y 2 + 2V we can write 

N-y 2 = 2f(P 2 ,I 2 ,y 2 ,L 2 ,yP) 

This equation implicitly defines y 2 as a function of P 2 , z 2 , N, 1? ,y ■ P, 
if we leave apart a very exceptional case where V linearly depends on 
y 2 in a special manner (such case is not realistic anyway). 
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Therefore it is sufficient that z 2 = const, for having also y 2 = const. , 
which in turn, according to (|69p . implies z ■ y = const.. 
Finally the five dynamical variables (|24p . or equivalently (|25p. remain 
constant on the circular orbit, [j 

It follows that A is a linear function of T on circular orbits. 

dV 

Theorem 1 If the interaction is such that \ z, V\ = and -7— — = 

d(z-y) 

there exist circular orbits; on these orbits the relative motion is periodic 
in terms of the center- of -mass time. 

dz 

Proof. According to (1441) we have that — = y and acording to (1451) 

d\ 

dy dV 

we get — = — 2-— - z. Equations (f44l(|4"5l) are similar to that of a 
dX oz 2 

three-dimensional two-body problem. The case when {z a , V} is zero 
corresponds to the classical problem of motion under a central force, 
where circular orbits are known to exist. 

As seen above, y 2 is constant. In the analogy between our system and 
that of Galilean mechanics, A plays the role of time and — y 2 represents 
the squared velocity. As well as a circular motion with a velocity of 
constant lenght is necessarily periodic in time, here we have that z and 
y are periodic functions of A. Since we consider circular motion, T is a 
linear function of A and vice versa, so periodicity in A implies periodicity 
inT. [] 

Example: any V(P 2 ,z 2 ) admits circular orbits. 



5 Extreme mass ratio, one-body limit 

We keep considering unipotential models. The case where one mass 
can be neglected in front of the other one is of practical interest when 
one tries to justify a resonable expression of V. Indeed it is naturally 
expected that in the limit of an extreme mass ratio we recover a system 
made of particle 1 moving in the external field created by particle 2, the 
latter undergoing rectilinear uniform motion. 

Without loss of generality we assume m\ <C TO2- Note that we cannot 
just put m2 to infinity. This can be seen already in the framework 
of Newtonian mechanics, because the gravitational potential created 
around particle 2 could not remain finite when 7712 tends to infinity. 
Therefore we shall rather put 



TTli = 7 7772 
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and study the limit for 7 — > 0. 

In order to alleviate calculations, let us set e = 7 2 . We have 

V = \ml{e + l) (70) 



v = \m 2 {e-l) (71) 



whence we derive 



» 2 -v 2 = \ml [(e + l) 2 — (e — l) 2 ] 

^ 2 - z/ 2 = e m 4 ( 72 ) 

We whish to investigate whether, looking at things from the rest frame, 
the center of mass E and the spacetime position %2 of the most heavy 
body actually coincide in the limit 7 — > 0. By formula (|20p we may 
write 

S = G + ( ? ^£)*-(^£)l/ (73) 

where Q = q2 + \z = q\ — \z. On the mass shell we have that P 2 = M 2 
so 

3 = *> + <k + w )z -i£ v (74) 

We can write m 2 = e m. 2 , , where, of course e = j 2 . 

We are interested in what happens when 7—7-0. 

In order to consider the most general case, let us set 

A = a ml (75) 

without assuming for the moment any restriction about the magnitude @ 
of a. 

In view of £70]) and ([73} we get 

(/i + A) 2 -v 2 = m\a (A + 2/x) + em\ 

(p + A) 2 — v 2 = m\ a [m\a + m\{e + 1)] + em\ 

(li + K) 2 - v 2 = m\ (a + l)(a + e) (76) 

Notice that 2(fi + A) = 77i|(l + e + 2a). Inserting into ([38]) yields the 
rigorous formula 

M 2 = m|[l + 2a + e + 2^(1 + a)(a + e)] (77) 

valid irrespective of the order of magnitude of a and e. 
Now we are in a position to make the following statement 



^ere we change the notation of [14] by suppression of a factor i. 
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Theorem 2 Provided we can neglect v}A| in front of m-i , we have 
that M 2 —7- — 2v, which entails that, at equal times, S and X2 coincide 
in the limit 7—^0. 



Indeed neglecting — ■ — - amounts to cancel a in (1771) . that yields — I as 

the limit of the ratio v/M 2 , making the second term in the right-hand 
side of (|74[) to vanish; remember that the third term vanishes at equal 
times, [j 

Owing to Propo. 3 the condition for this result is always satisfied for 
negative A. 

In contradistinction large positive values of A forbid E to coincide with 
the heavy body in the limit of an extreme mass ratio. In this case a > 
and formula (1771) can be written 



M 2 = m 2 2 [1 + 2a + e + 2^/(1 + a)a ^Jl + e/a] 

Define 

P = 2a + 2\/a 2 + a (78) 
Since a > it is clear that /3 > 2a. We get 

M 2 = m|(l + /3) + O(e) (79) 

now using ([7T|) yields 

\ + w> = 2ih) +0 « < 80 > 

which in general cannot vanish when e — > 0. [j 

This situation can be physically interpreted as follows: Considered at 
equal times, our covariant definition of the center of mass H reduces 
to that of Fokker and Pryce [16]. See also Moeller [18]. Accordingly 
we notice that H is in fact a center of energy; therefore not only the 
masses but also the energies must be taken into account. Even if mi is 
very small, it must be understood that we cannot consider particle 1 as 
a test particle when its motion involves a too large amount of energy. 



6 Toy Model 

. Consider a harmonic potential 

V = X "/P^z 2 (81) 
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with x a positive string constant (this potential differs from the one 
considered in [13] by because it is P 2 -dependent, which alows for the 
correct dimension of the coupling constant). 

The structure of the calculations derived from (|8ip is that of a nonrela- 
tivistic problem. Our notation is such that the relativistic potential and 
its nonrelativistic conterpart have opposite signs, so we have a positive 
A. 

We must compute F, as defined in (f52"jh according to ([5T]) and after 
solving the reduced equations of motion (|44 p ([35 ]) . The solution to this 
system is 



where A, B are mutually orthogonal spacelike constant vectors (they 
span the orbital plane, their lenghts are the half-axes of an ellipse) and 
C is a scalar constant; moreover we have 



which is always negative. It is clear that F will depend on A only 
through z 2 . Taking (|82p into account and fixing P a = k a we are left 
with 



z = A sin(OA + C) + B cos(OA + C) 
y = AQ cos(OA + C)- Bn sin(OA + C) 



(82) 
(83) 



Note that 




whence we derive 



< N >= 2 X M{A 2 + B 2 ) 



2 X M(a 2 + b 2 




setting A 2 



= —a 




b 2 . We obtain 



F = — 



X M[a 2 sm 2 (Q\ + C) + b 2 cos 2 (QX + C)] 



(84) 



In order to calculate / FdX we notice the primitive 



2 



/ (a 2 sin 2 (OA + C) + (b 2 cos 2 (OA + C)) dX = 

L 52 b 2 _ a 2 

A H — sin(OA + C) cos(OA + C) + const. 
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Finally we find 



r x a 2 + b 2 b 2 - a 2 

/ FdX = -xM [ — ±— A + 4 ^ sin(2f}A + 2C)] + const. (85) 

so there is a secular (linear) term plus a periodic correction; in this 
particular example r is a periodic function of A. 

Since z ■ P has a vanishing bracket with P 2 and z 2 it is obvious that 
{z ■ P, V} hence G, vanishes. 

M u 2 



Owing to (|34J) we are sure that — — > 0. But the question is 
dT 

about — . It can be directly read off (j84"|) that 
dX 

\F\ < xM(a 2 + b 2 ) 

r dT 

so the condition for -— > is 
dX 

4 ~ W > 1 1 

But 

\F\ < xM(a 2 + b 2 ) 
hence a sufficient condition 

AP _u^_ A 
4 ~ M 2> 2 



7 Conclusion 

Center of mass and relative motion are well understood concepts for iso- 
lated two-body systems, provided interaction is of the unipotential type 
and physical positions are fixed by the equal-time prescription. A large 
part of our picture can be made abstractly but most physical features 
show up in the light of the equal-time description. The main scheme 
was put forward many years ago [7], but several important consequences 
are considered here for the first time. 

In the present work we regard seriously the fact that the collective evo- 
lution parameter which arises in the reduced equations of motion (I44D 
(|45p is generally not a linear function of the center-of-mass time. This 
peculiarity (which doesnot concern the geometry of the orbit) automat- 
ically affects the schedule of relative motion. 

This point led us to consider an important exception, the case of circular 
orbits, where the relative motion is periodic in T. 
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Although the individual evolution parameters t\ , t 2 are generally not 
affine on the world lines, our equal-time treatment offers the possibility 
to end up with a description in terms of the center-of-mass time. 

On the other hand, examinating the case of an extreme mass ratio 
provides an illustration of the nature of center of mass in relativity. 
Indeed the result expressed in Theorem 2 forces one to interprete 5 as 
a center of energy rather than of mass; this is in agreement with an 
ancient literature and with the spirit of relativity. 

In this paper we considered relative motion essentially by analogy 
with a nonrelativistic one-body problem, naturally suggested by (I44D 
(|45p . But of course the question as to know whether (and how) a 
ficticious relativistic one-body system can be invoked, is relevant and 
deserves a separate publication. 



APPENDIX 

I. Reverting from fj,, v to mi,m2, formula (I38p implies this useful ap- 
proximation 



At first sight the appearence of tuq in this formula seems to indicate that 
the non-relativistic expression for reduced mas goes over to the relativis- 
tic realm without modification; notice however that we could replace mo 
by any positive m in this formula, provided that m = mo + 0(1/c 2 ) . 
This remarks leaves open the possibility that the "good" relativistic 
generalization of the reduced mass may coincide with mo only in the 

nonrelativistic limit, as happens for instance with Todorov's |19j reduced 

mim-2 

mass mt = . 

M 

II. Defining M a = P ■ p a /Mc 2 we have 




(86) 



But 2v = (m\ -m\)(? thus M x - M 2 

1 ; 

According to (f86|) we have — = 



Mi +M 2 = M, 



Mi - M 2 




2yP_2u 
Mc 2 ~ M 

2 



m\ + m-2 



1 



M ' 

+ 0(l/c 2 ), so finally 




+ 0(l/c 2 ) = mi - m 2 + 0(l/c 2 ) 



Finally M a = m a + 0(l/c 2 ). 
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